An improved exchange-correlation potential for polarizability and dissociation in DFT 
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We propose a novel approach to the problem of polarizabilities and dissociation in electric fields 
from the static limit of the Vignale-Kohn (VK) functional. We consider the response to the purely 
scalar part of the VK response potential. This potential has ground-state properties that notably 
improve over the full VK response density and over usual (semi-)local functionals. The correct 
qualitative behavior of our potentials means that it is expected to work well for polarizabilities in 
cases such as the H2 chain, and it will also correctly dissociate open-shell fragments in a field. 



Density functional theory (DFT) 0, H 0] has become 
the most popular electronic structure method in a wide 
range of problems in quantum chemistry achieving an 
unprecedented balance between accuracy and efficiency. 
But in some applications there are fundamental problems 
with standard exchange-correlation (xc) approximations 
such as LDA and GGAs. The response to electric fields is 
severely overestimated for long-range molecules and for 
molecular chains 0, [j| . Even for the simple H2 molecule 
dissociating in an electric field, LDA and GGAs yield un- 
physical fractional charges on each atom This is a 

problem in many topical applications, including molecu- 
lar electronics @, Q and nonlinear optics devices 0, [l(| ■ 

It is now well known that the functionals in these cases 
must depend in an ultranonlocal way on the density 
[Lllllli- There have been two approaches: The first is ex- 
act exchange (EXX): this has implicit non-local density- 
dependence through orbital-dependence [l3[. The sec- 
ond is Vignale-Kohn theory (VK) [TJ], where a linear 
response calculation is performed within time-dependent 
current-DFT (TDCDFT), a natural extension of time- 
dependent-DFT [l5| . VK uses functionals of the current- 
density: the (time-integrated) current through a small 
volume in space contains information about the density 
response far from that volume, so functionals that are 
local in the current-density are ultranonlocal in the den- 
sity. Both EXX and VK have field-counteracting terms 
that decrease the response compared with LDA/GGA, 
and so yield improved static polarizabilities for many 
polymer chains. An exception is the hydrogen chain, 
for which VK performs almost as poorly as LDA pd| . 
Another problem that needs non-local functionals is the 
dissociation of simple diatomics in electric fields. As the 
molecule is pulled apart, a field-counteracting step de- 
velops in the exact Kohn-Sham (KS) potential midway 
between the atoms. In the limit of large separation, the 
step size approaches a constant that realigns the highest 
occupied molecular orbitals of the two atoms in the case 
of two open-shell fragments, and vanishes in the case of 
two closed-shell fragments. EXX captures the step in the 
latter, but until now, no density functional approxima- 
tion has captured the step in the former case. VK has 



never before been applied to this problem. 

The calculations in utilise the zero-frequency limit 
of VK, which contains additional "dynamical" xc fields, 
on top of LDA [l|| . The VK xc vector potential has both 
longitudinal as well as transverse components. Static re- 
sponse calculations are technically outside the realm of 
validity of VK [lj], but its success in many such cases 
suggest such dynamical terms are also present in the true 
static functional and contain essential physics [ll], [ItJ • 

The fact that a transverse xc field persists in the 
static limit of VK means that, in contrast to usual re- 
sponse methods, caution must be used when interpret- 
ing the zero-frequency limit as a ground-state pertur- 
bation. Satisfaction of the adiabatic theorem requires 
that the response in the static limit is representable by 
a scalar KS potential [2(J: For a perturbation turned 
on slowly enough, the system remains arbitrarily close to 
the instantaneous ground-state. Functionals in any time- 
dependent theory should reduce to ground-state DFT 
ones, where xc effects are contained in a scalar potential. 
However, a transverse field, with its non-zero curl, repre- 
sents a non-conservative force, and so cannot correspond 
to a scalar potential underlying a conserved energy. 

In this paper we reconsider the VK response to electric 
fields from a new perspective. We first generate the self- 
consistent xc vector potential from a full VK calculation 
in the static limit, but then discard its transverse com- 
ponent. Thereby we eliminate the non-conservative part 
of the force. We gauge-transform the longitudinal part 
into a scalar potential which we view as a ground-state xc 
response potential. This approach is quite distinct from 
the previous use of VK , where the key player is the 
density response of the VK potential, n\ K . This is not 
the ground-state density response of the scalar poten- 
tial above, because n\ K is the full response to both the 
longitudinal and transverse fields of VK. Here, we con- 
sider the true ground-state response to just the scalar 
part of the VK potential. We show that this ground- 
state VK potential has desirable features arising from 
global field-counteracting terms. The dissociation limit 
of the electron-pair bond is correctly obtained, in con- 
trast to the notorious fractional charges that result from 
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all previous density functional approximations, including 

LDA, GGA and EXX HE El- F ° r CaS6S wher6 the VK 
response density has not performed well, eg. hydrogen 
chains, features of this scalar VK potential suggest that 
our approach will perform very well. Dynamical terms 
that proved crucial in the usual VK approach to polariz- 
abilities, play an even more fundamental role here. 

Consider a system initially in its field-free ground- 
state, of density Uo(r)- In VK response theory, xc- 
contributions to a perturbative time-dependent field are 
contained in a vector potential that is a local functional 
of the induced current-density j(ri) [lij . a xc [j](ri): in 
the frequency domain, 

z^a xc [j](r,c) = V<4 LDA (r) - V • *?(r, u)/n (r) (1) 



where, in the static limit, oj — > 0, 



Aioj 



drj dr{ 
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with u(r) = j(r)/no(r). Here /xc n ( r ) is a dynamical cor- 
rection to the static scalar xc kernel of the homogeneous 
electron gas /j™ in the long- wavelength limit, evaluated 
at the field- free density at r, i.e. 



lim/ s h r(«o(r)) 
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(3) 
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The static xc response of the homogeneous electron gas 
is given by the second-derivative of the xc energy density 
ex° m - This differs from the zero frequency limit of the 
scalar kernel by /xc" — 4/i X c n /3fi 2 , where /x X c n is the 
elastic shear modulus at oj — > jl9j. Thus VK contains 
dynamical corrections to ALDA that persist all the way 
through oj — > 0. The response kernels are defined in VK 
by taking the wavevector q of the perturbing field to zero 
before ta king the static limit, and these limits do not 
commute Hi EE E3, El ■ 

We now define a ground-state potential from the static 
limit of the VK response: Taking oj — * 0, the longitudinal 
component of Eq. |T]) is gauge-transformed to a scalar po- 
tential, Vxc [n] (r) , a non-local functional of the density: 
ioja. xc ^[j] (roj) = Vw X c [ n ] ( r ) • We restrict now to cylindri- 
cally symmetric linear systems, and approximate this po- 
tential in the following way. We consider Eqs. (0 and ([2]) 
for a purely one-dimensional inhomogeneity [1, El , so we 
effectively average over the transverse directions in our 
cases. We then obtain an approximate potential E El 
for linear systems: 



n'o(z) 



v&(z) = f*° m (z)Mz) - ^Kix d rW / dz'm(z') 

na{z) 

n' (7') 
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Here n' (z) — dno(r)/dz (where r = (0,0, z) is along the 
bond-axis) and ni(r) = V ■ j/(ioj) is the system's den- 
sity response, which is taken in the zero-frequency limit. 
This potential consists of four terms: the first is the LDA 
response, local in the density response. The second term 
is directly proportional to the local current-density re- 
sponse. The third and fourth display global behavior 
across the molecule, and are the key terms for the pur- 
poses of this paper. It follows from the structure of the 
fourth term that any polarization of the density, be it lo- 
cal or of charge-transfer nature, yields field-counteracting 
behavior. The third term tends to align along the field, 
but is generally smaller than the fourth. 

The first step is to run a zero-frequency VK response 
calculation on the chosen system in a weak external field, 
Ez, placed along the bond-axis, z. We utilize the ADF 
program package [2l| with the TDCDFT extension; sec 
Refs. [ll[ for implementation details. In the second step, 
the resulting density response is inserted in Eq. ([4| to de- 
fine our ground-state potential. We study three classes of 
systems that are challenging for usual semi-local density 
functionals: (i) dissociation of the electron-pair bond (ii) 
a dimer composed of two closed-shell fragments at large 
but finite separation, and (iii) a molecular chain. Atomic 
units are used throughout. 

(i) Dissociation of the electron-pair bond: H2-like sys- 
tems The step that forms in the exact KS potential when 
a molecule composed of open-shell units dissociates in an 
electric field prevents dissociation to fractionally charged 
species, and is consistent with the physical picture of two 
locally polarized species [5( . An analogous step occurs in 
field-free dissociation of a heteroatomic molecule com- 
posed of open-shell units E HI HI, Hi, HI . Its origin is 
static correlation, and it is particularly difficult for ap- 
proximations to capture, eluding not only LDA/GGAs 
but also EXX. Ref. [II has a density-matrix solution for 
this. 

The lower left panel of Figure [1] shows the field- 
free density (scaled by 0.01), the exact density response 
within bound-state perturbation theory, and the VK den- 
sity response n^ K for H2 at bond-length 10 au in a field of 
0.001 au. The exact response demonstrates local polar- 
ization with no charge-transfer, as expected. However, 
this is not true of the VK (or LDA) response density, 
which yields fractionally charged atoms. Correspond- 
ingly, the VK polarizability is grossly overestimated, as 
it is in the LDA. The top panel of this same figure shows 
the exact xc response potential (within bound-state 
perturbation theory) and the VK and LDA xc response 
potentials (i.e. subtracting the field- free potential). The 
exact was obtained by numerically inverting the KS equa- 
tion, for the exact KS bonding-orbital composed of po- 
larized atomic orbitals. The salient feature of the exact 
Uxc is the field-counteracting step which compensates for 
the difference in potential at the two separated atoms, 
and thus re-aligns them when the total potential is con- 
sidered. This step is missed by the LDA which has an 
along-field component, strengthening the applied field. 
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FIG. 1: H 2 at bond length 10 a.u., in a field of 0.001 a.u. 
The circles on the z-axis in the top panel represent the po- 
sition of the atoms, as in all subsequent figures. A similar 
field-compensating behavior was observed in other diatomics 
composed of open-shell fragments, eg. Li2. 



The step is also missing in EXX (not shown), as the po- 
tential there is simply minus half the Hartree. The VK 
potential does capture the step. In the right panel we see 
that the fourth term in Eq. ^ is responsible for this. The 
VK potential falls a little short of the exact step, but once 
the Hartree response potential (field-counteracting 0]) is 
added to the VK xc potential, we expect that the to- 
tal step completely compensates the change in potential 
created by the external field. 

The VK potential re-aligns the two atoms, and so the 
ground-state of this system would certainly not involve 
any charge transfer across the system, like in the exact 
case, in contrast to the VK density response (see Fig. 1). 
Indeed, as the bond-length R increases, one finds that 
the VK density globally polarizes even more, while the 
VK step size increases as ER, maintaining the atoms at 
the same level. As the molecule dissociates, stretched H2 
has the metallic-like feature that its HOMO-LUMO gap 
vanishes. This may underlie the reason why VK exactly 
captures the step, since VK is based on the response of a 
metallic system, the weakly inhomogeneous electron gas. 

It is important that the VK potential be evaluated 
on the "wrong" density response to the full VK fields: 
if evaluated on the exact response density, the field- 
counteracting term is reduced by a factor of about 20. 

(ii) Dimer of two closed-shell units Field-counteracting 
behavior also arises from the exchange interaction be- 
tween closed shells, as explained in Refs. 0, [26|. In con- 
trast to the step of (i), EXX methods can retrieve this 
step [1]. In Figure [2] we see that the step is also nicely 
reproduced by the VK potential, shown in the top panel. 
Also shown there is the exact xc response potential, cal- 
culated from highly accurate wavefunction techniques in 
Ref. [5[ . Although details of the exact potential are miss- 
ing in the VK response potential, the step is clearly cap- 
tured, with the correct magnitude of drop in the potential 



FIG. 2: The H2-dimer, with H2- H2 distance 3 a.u., and 
"intramolecular" H-H bondlength 1.4 a.u, in a 0.001 au field. 
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FIG. 3: H 2 chain in a field of 0.005au: OEP (solid), KLI 
(dotted), LDA (dashed), and VK potentials (+). The H-H 
distances alternate between 2 and 3 a.u. 

between the up- and down-field molecules. In contrast, 
the LDA or any GGA has no net drop. 

The right panel shows that the net step in the VK 
potential is a result of competition between the third and 
fourth terms. Note that in this case there is largely local 
polarization in the molecules, with small charge transfer 
(hardly discernable in the lower left panel). 

The VK density response and polarizability are again 
very close to the LDA ones, both overestimating the ex- 
act polarizability. However, the ground-state response 
using the VK potential shown, is expected to yield po- 
larizabilities much closer to the exact one, because of the 
field-counteracting behavior. 

The magnitude of the exact step between two closed- 
shell units decreases as the intermolecular separation in- 
creases [H . This is not the case for the VK step asymptot- 
ically. Essentially, the VK step arises from integrals over 
the current-density (fourth term in Eq. ([!])); this integral 
persists at large separation, due to the local polarization. 

(Hi) Long-chain molecules Polarizabilities and hyper- 
polarizabilities of long-chain polymers are notoriously 
overestimated in LDA/GGA calculations. Recent work 
has shed much light on the role of ultranonlocal density 
dependence in these systems. It has been demonstrated 
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that orbital functionals such as EXX within OEP can 
significantly improve the polarizabilities 0, [^, H3 |. due 
to field-counteracting behavior in the exchange poten- 
tial. Approximate exchange methods such as KLI (2?| 
also give corrections over LDA, although not as strongly 
as full exchange @. For example, for the hydrogen chain 
in Figure [21 the LDA, KLI, and OEP polarizabilities are 
114.6, 90.6, and 84.2 a.u, respectively. The VK func- 
tional was also applied [ll| and gave 110.20 a.u, hardly 
an improvement over LDA; in contrast to its dramatic 
improvement for many other molecular chains. Again, 
the scalar VK potential has field-counteracting behavior, 
as shown in Figure [3) This potential will give improved 
results (comparable to KLI) when used in a computation 
of ground-state response and polarizabilities. 

In Ref. 0] , the importance of intermolecular barriers in 
the field-free ground-state potential was stressed: these 
are lacking in LDA, underestimated in KLI, but captured 
well in OEP. Although formal arguments point to LDA as 
the correct field-free potential to be used in conjunction 
with the VK response potential, it will be interesting to 
compare polarizabilities using the VK response potential 
on top of a field-free OEP potential. 

In summary, we have shown that the scalar part 
of the VK response potential contains crucial field- 
counteracting terms making it a good candidate for re- 
sponse properties of effectively one-dimensional systems 
in electric fields. It is the first approximation that cap- 
tures the step in a molecule composed of open-shell frag- 
ments dissociating in an electric field. It captures field- 
counteracting terms in systems of two or more closed- 
shell units, and is promising for long-chains, being nu- 



merically less intensive than EXX methods which also 
have been successful for these problems. The scalar VK 
potential works well even when the full VK response does 
not! In all cases studied, n\ K is very close to the LDA 
response, suggesting that the effects of the transverse VK 
and dynamical-longitudinal fields somewhat cancel. Sat- 
isfaction of the adiabatic theorem means that the trans- 
verse component of the response field should vanish in the 
static limit, suggesting that VK should be corrected by 
dropping its transverse part. This supports our results, 
where the ground-state response to the longitudinal com- 
ponent was considered. An alternative choice is to dis- 
card the transverse field from the start, and consider the 
self-consistent response to the purely longitudinal field. 
The results are not as good: the field-counteracting na- 
ture of Eq. (fj| is most effective when evaluated on the 
full (wrong) density-response n\ K rather than that ob- 
tained self-consistently. The non-self-consistent aspect 
of our present approach is somewhat dissatisfying from 
a rigorous viewpoint; on the other hand, there are many 
other situations in DFT where such post-self-consistent 
approaches have success eg. Ref. (28|. 

Finally, we note that the VK response potential smears 
over some local details of the true response potential (see 
e.g. Fig. [2]). Work is underway to investigate if this has 
a significant effect on the global polarizability. 
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and National Science Foundation Career CHE-0547913. 



[1] P. Hohenberg and W. Kohn, Phys. Rev. 136, B 864 
(1964). 

[2] W. Kohn and L.J. Sham, Phys. Rev. 140, A 1133 (1965). 
[3] W. Kohn, Rev. Mod. Phys. 71, 1253 (1999). 
[4] B. Champagne et. al, J. Chem. Phys. 109, 10489 (1998). 
[5] S.J. A. van Gisbergen et. al., Phys. Rev. Lett. 83, 694 

(1999); Phys. Rev. A 62 012507 (2000). 
[6] O.V. Gritsenko and E.J. Baerends, Phys.Rev. A 54, 1957 

(1996). 

[7] C. Toher et al., Phys. Rev. Lett. 95, 146402 (2005). 

[8] N. Sai et al., Phys. Rev. Lett. 94, 186810 (2005). An effec- 
tively one-dimensional approach was also adopted here. 

[9] S. Kummel, L. Kronik, and J. P. Perdew, Phys. Rev. 
Lett. 93, 213002 (2004). 
[10] D.R. Kanis, M.A. Ratner, and T.J. Marks, Chem. Rev. 
94, 195 (1994). 

[11] M. van Faassen et al., Phys. Rev. Lett. 88, 186401 (2002); 

J. Chem. Phys. 118, 1044 (2003). 
[12] P. Mori-Sanchez, Q. Wu, and W. Yang, J. Chem. Phys. 

119, 11001 (2003). 
[13] J.D. Talman and W.G. Shadwick, Phys. Rev. A 14, 36 

(1976). 

[14] G. Vignale and W. Kohn, Phys. Rev. Lett. 77, 2037 
(1996); G. Vignale, C.A. Ullrich, and S. Conti, Phys. 



Rev. Lett. 79, 4878 (1997). 
[15] E. Runge and E.K.U. Gross, Phys. Rev. Lett. 52, 997 

(1984) ; 

[16] C. A. Ullrich and G. Vignale, Phys. Rev. B, 58, 15756 

(1998); ibid 65, 245102 (2002). 
[17] G. F. Giuliani and G. Vignale, Quantum Theory of the 

Electron Liquid, Cambridge, 2005. 
[18] E.J. Baerends, Phys.Rev. Lett. 87, 133004 (2001); M.A. 

Buijse and E.J. Baerends, Mol. Phys. 100, 401 (2002). 
[19] Z. Qian and G. Vignale, Phys. Rev. B. 65, 235121 (2002). 
[20] E.K.U. Gross, J.F. Dobson, and M. Petersilka, in Topics 

in Current Chemistry, 181, 81 (1996). 
[21] ADF 2003.01, available from |http: / /www.scm.com| 
[22] J.P. Perdew et al, Phys. Rev. Lett. 49, 1691 (1982). 
[23] J.P. Perdew and M. Levy, Phys. Rev. B 56, 16021 (1997). 
[24] J.P. Perdew, in Density Functional Methods in Physics, 

ed. R.M. Dreizler and J. da Providencia (Plenum, NY, 

1985), p. 265. 

[25] CO. Almbladh and U.von Barth, Phys. Rev. B. 31, 3231, 

(1985) . 

[26] O.V. Gritsenko and E.J. Baerends, Phys. Rev. A 64, 
042506 (2001). 

[27] J.B. Krieger, Y. Li, and G. J. Iafrate, Phys. Rev. A 46, 
5453 (1992). 



[28] J. P. Perdew et al., Phys. Rev. Lett. 82 2544 (1999). 



